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1.  Introduction 


“Yes,  we  have  to  divide  up  our  time  like  that,  between 
our  politics  and  our  equations.  But  to  me  our  equations  are 
far  more  important,  for  politics  are  only  a  matter  of  present 
concern.  A  mathematical  equation  stands  forever." 

Attributed  to  Albert  Einstein 


This  is  a  continuation  of  the  previous  reports  in  this  series  ( 1 ,  2),  where  the  motivation  is 
presented.  In  those  reports,  we  gave  several  examples  of  generating  complete 
thermodynamically  consistent  equations  of  state  (EOS).  The  methodology  used  there  was  based 
on  combining  some  ad-hoc  generalizations  of  classical  models  with  available  experimental  data. 
The  EOS  have  been  presented  in  analytical  form  as  functions  of  the  experimental  data.  Their 
numerical  implementation  has  been  discussed  by  Bilyk  et  al.  (J). 

We  begin  with  the  general  analysis  which  does  not  rely  on  any  particular  assumptions  but  those 
of  classical  thermodynamics  (like,  for  instance,  the  assumption  of  dealing  with  a  two-parameter 
system).  In  this  respect,  our  method  reminds  the  remarkable  short  letters  of  Peek  and  Salsburg 
( 4 )  and  of  Zel’dovich  (5),  although  the  differences  between  the  two  methods  are  quite  significant 
and  will  be  discussed  elsewhere.  Unfortunately,  so  far  the  Peek-Salburg-Zel’dovich  approach 
did  not  find  many  followers  neither  in  the  former  Soviet  Union  nor  in  the  West.  The  only 
exception  that  this  author  is  aware  of  is  the  paper  of  Fortov  and  Krasnikov  (6). 

In  general,  our  approach  leads  to  a  pair  of  integro-differential  equations  which,  probably,  should 
be  solved  numerically.  We  then  proceed  with  the  special  case  of  constant  heat  capacity.  This 
special  case  is  particularly  important  for  applications  and,  fortunately,  permits  explicit  solution. 
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2.  Required  Notions  and  Basic  Formulas 


In  this  report,  we  present  a  methodology  which  does  not  rely  on  any  a  priori  guesses  regarding 
the  structure  of  the  EOS.  The  only  assumption  made  is  that  the  substances  under  study  can  be 
treated  as  two-parameter  thermodynamic  systems  with  sufficiently  smooth  EOS. 

According  to  our  definition,  for  a  medium  with  two  thermodynamic  parameters,  the  complete 
EOS  for  hydrocode  is  nothing  but  the  specific  entropy  density  S  presented  as  a  function  of  the 
specific  volume  V  internal  energy  density  E  : 

S=S(V,E)  (1) 


In  terms  of  this  function,  the  pressure  P=  P(V,  E  )  and  the  absolute  temperature  can  be 
presented  in  the  following  form: 

dS(V,E) 

P(V’E)  =  dS(y  ^  ’  T(V>E )  =  dS(y  Ej  • 

6E  dE 


The  isentrope  in  the  ( V,E )  parametric  space  is  a  curve  E  =  Es  (V )  along  which  entropy  density 
remains  constant.  We  use  notation  E  =  E*s*  (V )  for  the  isentrope  passing  through  the  state 
(V*,  E'  Y  In  figure  1,  two  isentropes,  E  =  E**  ( V )  and  E  =  Ef  ®  (V )  ,  are  presented  in  green:  the 
former  passing  through  the  state  (y\E^  and  the  latter  passing  through  the  state  (V  s,  E'  ;  j . 
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Given  the  starting  state  (v*,E* ) ,  the  Hugoniot  adiabat  E  =  E**  (V )  is  the  curve  presented  all  the 
states  ( V,E ) ,  which  can  be  reached  from  the  starting  state  by  means  of  the  shock  waves  of 
different  intensity.  In  the  starting  state  (v \E ),  the  isentrope  and  the  Hugoniot  adiabat  have 
common  tangent  lines.  However,  somewhat  farther  from  that  state,  the  Hugonuit  adiabat  has  a 
bigger  tangent  line  and,  quite  often,  it  even  has  a  vertical  asymptote  at  a  finite  value  of  the 
specific  volume.  There  are  two  more  principal  distinctions  between  isentropes  and  Hugoniot 
adiabats.  First,  each  isentrope  can  be  extended  from  its  pole  (v®  ,E®^  in  both  directions;  the 
Hugoniot  adiabat  can  be  extended  in  a  single  direction  only.  Second,  if  a  new  point  (v°,Ee^j 
on  the  Hugoniot  adiabat  E  =  E**  (V )  is  chosen  as  a  starting  point,  then  the  new  Hugoniot  adiabat 
E  -  E®°  (V )  does  not  coincide  with  E  =  En  (V ) .  Further  properties  of  the  Hugoniot  adiabat  can 
be  found  in  numerous  textbooks,  including  the  classical  ones  (7,  8). 

Given  the  isentrope  E  =  Ef®  (V ) ,  the  pressure  values  P  =  P®®  (V )  on  this  very  isentrope  can  be 
calculated  by  means  of  differentiation  only  as 


_  dE**(v ) 
K  ’  dV 


(3) 
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The  opposite  is  also  true.  Given  the  pressure  values  P  =  P®®  (V ) ,  the  equation  of  the  isentrope 
can  the  calculated  by  means  of  straightforward  integration: 

Er{v)=\dtpr{t).  (4) 

V 

Experimenters  usually  measure  the  pressure  P  =  P®®  (V )  on  isentropes,  not  the  energy 
density  E  =  E®®  (V )  .  When  dealing  with  the  Hugoniot  adiabat,  the  situation  is  opposite  (9, 10). 
The  first  things  experimenters  measure  in  typical  experiments  on  shock  waves  are  the  velocities 
of  the  material  particles  in  front  of  and  behind  the  shock  wave  and  the  velocity  of  the  shock  front 
itself.  Then,  using  those  measurements  and  the  conservation  equations  across  the  shock  front,  it 
is  possible  to  calculate  the  specific  energy  density  E  =  E**  (V )  and  the  pressure  P  =  PH  (V ) 
behind  the  shock  front.  The  functions  E  =  EH  (V )  and  P  =  P^  (V )  characterizing  the  Hugoniot 
adiabat  are  interrelated  according  to  the  classical  Rankine-Hugoniot  formula 


P?(V) 


2E”(V)-2E*  +P*  (V-V*) 
V-V* 


(5) 


Summarizing,  having  the  experimentally  determined  data  on  P  =  P®@  (V )  and  E  =  E"H  (V ) ,  we 
can  automatically  calculate  two  more  functions:  E  =  E®®  (V )  and  P  =  P"  (V ) . 

If  we  have  the  experimentally  determined  specific  heat  capacity  Cv  =C(V,E )  at  fixed  volume, 
then  the  complete  EOS  for  hydrocode  can  be  presented  in  the  following  form: 


E 

s(y,£)  =  A(y)+J 

E* 


els 

H(V,^)  +  A(V)  ’ 


(6) 


where  the  function  H(V,£’)  is  defined  as 


E 

H  (v,E)=jdr,q;1 

E* 


(7) 


whereas  the  functions  A(V)  and  A(V)  should  be  recovered  from  additional  experimental  data. 


The  value  of  the  function  A(E)  at  V  =  V  is  equal  to  the  value  of  the  absolute  temperature  T 
in  the  state  (v*,E*y 

A(v*)=T*  .  (8) 
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3.  The  EOS  Recovery  for  a  General  Heat  Capacity  Function  Cv  (V,E) 


The  recovery  of  the  EOS  with  non-constant  heat  capacity  is  important  for  application  (see 
Segletes  [ 11 ]  and  references  therein).  In  order  to  recover  the  complete  EOS  (equation  6),  we 
have  to  express  the  functions  A(E )  and  A  (E )  in  terms  of  experimental  data.  To  find  those 

functions  from  the  isentrope  E  =  Ef®  (V )  and  the  Hugoniot  adiabat  E  =  En  (V  )  actually  means 
to  express  the  functions  A(E)  and  A  (E)  in  terms  of  the  functions  Cv  =C (E,  E),  E  =  E®®  (E ) , 
and  E  =  E ^  ( E )  .  From  the  previous  discussion,  it  follows  that  we  automatically  get  two 
additional  functions,  P  =  E®®  (E )  and  P  =  Pn  ( E ) ,  which  can  be  calculated  with  the  help  of 
equations  3  and  5,  respectively. 


In  fact,  the  functions  A(E)and  A(E)  can  be  recovered  with  the  help  of  the  function 

E  =  E®®  (V )  and  E  =  EH  ( E )  by  means  of  solving  the  following  pair  of  the  integro-differential 

equations: 


A(V)+  | 


•oo 


H(E,£)  +  A(E) 


-S®  =  0 


(9) 


and 


_d_ 

clV 


Ef®(V) 


Eh{V) 


d% 


dE*H  (E) 


dV 


+  Ph\V) 


H(E,^)  +  A(E)  H(E,E;(E))  +  A(E) 


=  0 


(10) 


that  should  be  combined  with  the  initial  condition  (equation  8). 

In  fact,  the  system  (equations  9  and  10)  can  be  solved  sequentially.  First,  we  solve  equation  10 
for  the  function  A(E ) ,  and  then  we  calculate  the  function  A(E )  with  the  help  of  equation  9. 

The  solutions  should  be  then  inserted  into  equation  6. 
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4.  Explicit  Recovery  Formulas  for  the  Case  of  Constant  Heat  Capacity 


Solution  of  the  general  system  (equations  9  and  10)  requires  computer-based  methods. 
Fortunately,  an  explicit  solution  is  possible  in  the  case  of  constant  heat  capacity.  This  special 
case  is  important  in  many  practical  applications.  In  addition,  when  implementing  computer- 
based  methods,  explicit  solutions  play  a  crucial  role  in  validation  and  verification  procedures. 

In  this  case,  the  solution  of  equation  10  with  the  initial  data  (equation  8)  for  function  A  ( V7 )  can 
be  presented  in  the  following  explicit  form: 

A(y)  =  r*[i-^(y,v*)].  (ii) 

Here,  and  in  the  following,  we  need  two  non-dimensional  functions  e and  r(l/,V  j of 
two  variables  which  can  be  presented  in  the  form 

V 

1  ^  — j* 

£(v,V*)  =  —  \d07r(0)T(V,O),  z(V,6)  =  e°  ,  (12) 

in  which  two  non-dimensional  functions  uj{V),  t(V)  of  one  variable  are  expressed  in  terms  of 
the  measurable  quantities  only  as 

nlv)  v-P»(v)-rr(v)  V  P?(v)Er(v)-rr(v)EZ(v) 


e;*(V)-£“(V) 


E~(V)-Ef®{V) 


In  the  case  of  constant  heat  capacity,  the  skeleton  (6)  reduces  to  the  simpler  form 


s(v  F\-^  =  C  I„  E-E'  +  Cv\(v) 

'  ’  )  v  £f’(V)-£‘  +  C„A(V)  ' 


Inserting  the  solution  (equation  11)  in  the  formula  (equation  14),  we  arrive  at  the  following 
complete  equation  of  state  for  hydrocode  (where  we  made  a  special  choice  E  =  CT  for  the 
arbitrary  constant) 


S(V,E)-S°  =Cv\n 


E-cvrs(vy ) 

E^{V)-CvT*s{vy )  ' 
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Associated  with  equation  15,  expressions  of  the  absolute  temperature  and  pressure  are  as 
follows: 


T(V,E)  = 


E-CvT*e(v,V* ) 


G, 


(16) 


and 


P  - 

(V)  =  H  s 


F  —  f 


rpQ9Q9  p**  _  E1  pQ9 

F  —^S  rH  rS 


F  —  F 


(17) 


The  formula  of  the  pressure  (equation  17)  shows  that  P(V,E )  is  a  linear  function  of  the  internal 

energy  density,  as  it  should  be  in  the  case  of  constant  heat  capacity.  In  fact,  the  formula  (17)  is 
more  general,  and  it  does  not  require  the  assumption  of  constant  heat  capacity.  It  is  valid  when 
the  heat  capacity  is  arbitrary  function  of  entropy.  In  many  applications  of  the  hydrocode  in 
ballistics,  only  the  incomplete  EOS  (equation  17)  is  required,  not  the  complete  EOS  (equation 
15). 


5.  Conclusion 


We  analyzed  the  important  problem  of  thermodynamically  consistent  recovery  of  the  complete 
EOS  for  hydrocode  from  experimental  data.  More  specifically,  we  demonstrated  how  to 
accomplish  the  recovery  based  on  the  measurements  of  the  heat  capacity  Cv  =C  (V,  E)  at 
constant  volume,  of  the  isentrope  Es  =  £®®  (V ) ,  passing  through  the  state  (V®,E®j,  and 
Hugoniot  adiabat  EH  =  E**  (V )  with  the  initial  state  (  V  \E  j . 

We  do  not  make  any  special  physical  hypothesis  about  the  EOS  to  be  recovered. 

We  assume  that  our  media  can  be  adequately  described  as  a  two-parameter  model,  that  the 
generated  EOS  satisfy  the  thermodynamic  stability  conditions,  and  that  phase  transformations 
can  be  ignored  (see  Grinfeld  [72]  for  further  discussion). 

It  is  demonstrated  that  the  recovery  basically  reduces  to  solving  a  nonlinear  integro-differential 
(equation  10).  Generally  speaking,  this  equation  cannot  be  solved  analytically.  At  the  same 
time,  it  allows  quite  a  straight-forward  implementation  into  a  computer-based  numerical 
analysis. 

In  the  practically  important  case  of  constant  heat  capacity,  the  integro-differential  equation  10 
allows  explicit  solution.  The  solution  is  given  by  the  formulas  in  equations  11  and  12.  The 
corresponding  complete  EOS  for  hydrocode  is  given  by  the  explicit  formula  (equation  15).  The 
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incomplete  EOS  for  the  absolute  temperature  and  pressure  are  given  by  the  formulas  in  equations 
16  and  17,  respectively. 

The  procedure  of  the  recovery  of  the  EOS  from  the  measured  constant  heat  capacity  C  at 
constant  volume,  measured  pressure  P00  (V )  on  the  isentrope,  and  the  internal  energy  density  on 
the  Hugoniot  adiabat  E**  (V )  is  the  following.  First,  using  the  functions  P00  (V )  and  E**  (V ) 
and  the  relationships  (equations  4  and  5),  we  calculated  two  additional  functions  Ef 0  (V )  and 
Pn  (V )  .  Then,  using  the  definitions  (equation  13),  we  calculate  the  functions  kt(E)  and  7i(V) . 
Afterwards,  using  the  definitions  (12),  we  calculate  the  functions  s(v, V”)  and  ^  .  At 

last,  we  substitute  the  function  s{V  ,V  )  into  the  relationships  (equations  15-17)  in  order  to 
determine  the  complete  EOS  for  the  entropy  potential  S  (V,E)  and  two  incomplete  EOS  for  the 
absolute  temperature  T(V,E )  and  pressure  P(V,E) . 

An  important  issue  of  tabulating  the  presented  analytical  forms  for  the  EOS  is  discussed 
elsewhere  (13, 14). 
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